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Workpiece  material  flow  and  stirring/mixing  during  the  friction  stir  welding  (FSW)  process  are  investigated 
computationally.  Within  the  numerical  model  of  the  FSW  process,  the  FSW  tool  is  treated  as  a  Lagrangian 
component  while  the  workpiece  material  is  treated  as  an  Eulerian  component.  The  employed  coupled 
Eulerian/Lagrangian  computational  analysis  of  the  welding  process  was  of  a  two-way  thermo-mechanical 
character  (i.e.,  frictional-sliding/plastic-work  dissipation  is  taken  to  act  as  a  heat  source  in  the  thermal- 
energy  balance  equation)  while  temperature  is  allowed  to  affect  mechanical  aspects  of  the  model  through 
temperature-dependent  material  properties.  The  workpiece  material  (AA5059,  solid-solution  strengthened 
and  strain-hardened  aluminum  alloy)  is  represented  using  a  modified  version  of  the  classical  Johnson-Cook 
model  (within  which  the  strain-hardening  term  is  augmented  to  take  into  account  for  the  effect  of  dynamic 
recrystallization)  while  the  FSW  tool  material  (AISI  H13  tool  steel)  is  modeled  as  an  isotropic  linear-elastic 
material.  Within  the  analysis,  the  effects  of  some  of  the  FSW  key  process  parameters  are  investigated  (e.g., 
weld  pitch,  tool  tilt-angle,  and  the  tool  pin-size).  The  results  pertaining  to  the  material  flow  during  FSW  are 
compared  with  their  experimental  counterparts.  It  is  found  that,  for  the  most  part,  experimentally  observed 
material-flow  characteristics  are  reproduced  within  the  current  FSW-process  model. 


Keywords  friction  stir  welding,  material-flow  simulation  and 
analysis,  process  modeling 


1.  Introduction 

Friction  stir  welding  (FSW)  is  a  relatively  new  solid-state 
material  joining  process  which  was  invented  at  The  Welding 
Institute,  Cambridge,  UK  in  1991  by  Thomas  (Ref  1).  Within 
this  joining  process,  a  non-consumable  hard-material  welding 
tool  is  used  to  generate  (via  frictional-sliding  and  plastic-work 
dissipation)  sufficient  amount  of  heat,  in  the  workpiece  material 
surrounding  the  tool/workpiece  interface,  required  for  success¬ 
ful  welding.  A  detailed  description/analysis  of  the  FSW  process 
could  be  found  in  the  seminal  work  by  W.  M.  Thomas  and  co¬ 
workers  (Ref  1,  2)  and,  hence,  will  not  be  provided  in  this 
article. 

While  it  has  been  demonstrated  that  a  number  of  metallic 
and  polymeric  materials  can  be  joined  using  FSW,  the  industrial 
interest  has  been  primary  in  the  welding  of  aluminum  alloys. 
For  a  wide  variety  of  aluminum-alloy  grades  (even  for  those 
previously  considered  as  unweldable),  it  was  clearly  demon¬ 
strated  that  high-quality/defect-free/high-strength  FSW  joints 
(in  a  1-50  mm  thickness  range)  can  be  produced.  The 
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additional,  often-cited  advantage  of  FSW  is  its  ability  to  be 
carried  out  in  various  positions  (e.g.,  horizontal,  vertical, 
overhead,  orbital,  etc.).  In  regard  to  different  weld  designs, 
FSW  is  most  often  used  (but  is  not  limited)  in  the  production  of 
butt  and  lap  joints. 

It  has  been  generally  established  that  the  FSW-tool  design  is 
an  important  FSW  process  parameter  affecting  all  aspects  of  the 
joining  process  (e.g.,  temporal  evolution  and  spatial  distribution 
of  the  material  and  heat  transport/flow,  material  microstructure 
and  properties,  the  presence/absence  of  weld  flaws/defects, 
etc.).  Typically,  an  FSW-tool  consists  of  three  main  sections, 
Fig.  1  (a)  shank;  (b)  shoulder;  and  (c)  pin.  Furthermore,  the 
FSW  process  is  generally  carried  out  using  a  conventional 
milling  machine,  within  which  the  FSW-tool  is  rotated  about  its 
(longitudinal)  axis,  while  the  workpiece  is  held  in  place  using 
a  fixturing  device.  The  main  role  of  the  shank  is  to  provide 
a  connection  to  the  milling  machine  spindle  and,  thus,  to  apply 
a  torque  to  the  FSW-tool.  The  primary  role  of  the  tool  shoulder 
is  to  apply  the  required  level  of  the  downward  contact  pressure 
onto  the  workpiece,  facilitate  generation  of  the  heat  and  ensure 
efficient  stirring  of  the  material  and  to  confine  the  thermally 
softened  workpiece  material  and  prevent  its  escape  from  the 
weld  region.  To  ensure  the  desired  level  of  the  tool-shoulder 
functionality,  the  bottom  face  of  the  shoulder  generally 
possesses  an  upright  truncated-conical  profile  as  well  as 
scrolls/spirals.  As  far  as  the  pin  is  concerned,  its  main  role  is 
to  provide  an  additional  contribution  to  the  heat  generation  (via 
frictional-sliding  and  plastic-work  dissipation)  and  to  ensure 
efficient  and  thorough  mixing  of  the  workpiece  material  in  the 
longitudinal,  transverse,  and  vertical  directions.  To  improve  the 
functionality  of  the  FSW-tool,  the  pin  is  often  tapered  and 
contains  external  features  such  as  threads,  flutes,  flats,  etc. 
Since  the  extents  of  material  mixing  and  heating  below  the  tool- 
pin  are  quite  low  and  mainly  localized  in  a  thin  region  adjacent 
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to  the  pin  bottom,  the  length  of  the  pin  is  generally  chosen  in 
such  a  way  that  it  extends  to  a  depth  very  close  to  the 
workpiece  bottom  face.  This  ensures  that  the  FSW  weld  is 
complete  and  defect  free. 

A  schematic  of  a  typical  FSW  set-up  used  in  the  fabrication 
of  butt-joints  is  displayed  in  Fig.  2.  The  FSW  fabrication  of 
butt-joints  generally  involves  the  following  steps:  (a)  The 
workpiece  components  with  mating  surfaces  are  subjected  to 
high  contact  pressures  and  secured  to  a  rigid  backing  plate 
using  a  fixturing  device.  This  prevents  potential  separation, 
sliding,  or  lifting  of  the  workpiece  components  and  enables 
positioning  of  the  work-piece  material  at  an  angle  relative  to  the 
tool  longitudinal  axis;  (b)  The  welding  tool  secured  within  the 
milling  machine  tool  holder  is  spun  to  the  desired  rotational 
speed  and  is  slowly  plunged  into  the  workpiece  components 
contact-interface  region  until  the  tool  shoulder  contacts  (and 
slightly  indents)  the  work-piece  upper  surface  and  the  lowest 
point  of  the  pin  is  very  near  the  workpiece  back  face;  (c)  The 
rotating  tool  is  then  forcibly  advanced  along  the  workpiece 
components  contact  interface  while  the  workpiece  material  is 
stirred,  extruded  around  the  tool,  and  subsequently  forged  in  the 
wake  of  the  tool;  and  (d)  At  the  completion  of  the  welding 
process,  the  tool  advancement  is  terminated  and  the  tool  is 
retracted  (i.e.,  pulled  out  of  the  workpiece)  while,  the  tool 
continues  to  rotate.  Once  the  tool  is  fully  retracted,  the  spindle 
rotation  is  terminated  and  the  weldment  is  removed  from  the 


fixturing  device.  Clearly,  a  hole  (which  needs  to  be  eliminated) 
is  left  in  the  weldment  upon  retraction  of  the  FSW  tool. 

Since  the  FSW  tool  rotates  in  the  same  (clockwise  or 
counterclockwise)  direction  during  welding,  the  thermo¬ 
mechanical  and  microstructural  fields  associated  with  the 
FSW  process  are  not  symmetric.  Consequently,  the  side  of 
the  weld  on  which  the  tangential  component  of  the  rotational 
velocity  and  the  translational  velocity  of  the  FSW  tool  are  of 
the  same  sense  is  commonly  referred  to  as  the  advancing  side 
while  the  other  side  of  the  weld  is  called  the  retreating  side.  In 
addition,  the  following  terms  are  often  used  to  denote  various 
regions  of  the  weld  relative  to  the  position  of  the  tool:  (a) 
Regions  ahead  and  behind  the  tool  are  referred  to  as  the  leading 
and  trailing  sides,  respectively;  and  (b)  the  top  and  bottom  of 
the  weld  are  denoted  as  the  crown  and  root,  respectively. 

FSW  normally  involves  complex  interactions  and  competi¬ 
tion  between  various  mass  and  heat  transport  phenomena, 
plastic  defonnation  and  damage/fracture  mechanisms,  and 
microstructure  evolution  processes  (Ref  3-10).  Consequently, 
the  material  microstructure  (and  mechanical  properties)  in  the 
weld  region  are  highly  complex  and  spatially  diverse.  Metal- 
lographic  examinations  of  the  FSW  joints  typically  reveal  the 
existence  of  the  following  four  weld  zones  (not  counting  the 
unaffected/base-metal  zone  which  is  far  enough  from  the  weld 
so  that  material  microstructure/properties  are  not  altered  by  the 
joining  process)  (Fig.  3)  (a)  the  heat-affected  zone  (FIAZ)  in 
which  material  microstructure/properties  are  effected  only  by 
the  thermal  effects  associated  with  FSW.  While  this  zone  is 
normally  found  in  the  case  of  fusion-welds,  the  nature  of  the 


Trt«l  ShmiMor  U/iHfh 


Fig.  1  Typical  friction  stir  welding  (FSW)  tool  used  for  joining  of  Fig.  3  A  schematic  of  the  main  microstructural  zones  associated 
aluminum  alloys  with  the  typical  FSW  joint 


Fig.  2  A  schematic  of  the  friction  stir  welding  (FSW)  process 
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microstructural  changes  (e.g.,  recovery,  recrystallization,  grain 
growth,  precipitate  coarsening,  precipitate  dissolution,  etc., 
depends  on  the  workpiece  material  chemical  composition  and 
microstructural  state)  may  be  different  in  the  FSW  case  due  to 
generally  lower  temperatures  and  a  more  diffuse  heat  source; 
(b)  the  thermo-mechanically  affected  zone  (TMAZ)  is  located 
closer  than  the  HAZ  to  the  workpiece  components  contact 
interface.  Consequently,  both  the  thennal  and  the  mechanical 
aspects  of  the  FSW  process  affect  the  material  microstructure/ 
properties  in  this  zone.  While  the  original  grains  are  retained  in 
this  zone  they  experience  a  considerable  amount  of  plastic 
defonnation  during  the  FSW  process.  Consequently,  it  is  often 
observed  (using  metallographic  analysis  of  the  transverse  and 
horizontal  weld  sections)  that  the  grains  have  been  subjected  to 
bending;  (c)  the  weld  nugget  is  the  innermost  zone  of  an  FSW 
joint.  As  a  result  of  the  way  the  material  is  transported  from  the 
regions  ahead  of  the  tool  to  the  wake  regions  behind  the  tool, 
this  zone  typically  contains  the  so  called  onion-ring  features. 
The  material  in  this  region  has  been  subjected  to  most  severe 
conditions  of  plastic  deformation  and  high-temperature  expo¬ 
sure  and  consequently  contains  a  very-fine  dynamically  recrys¬ 
tallized  (equiaxed  grain  microstructure)  weld  microstructure. 
The  width  of  the  nugget  is  slightly  larger  than  the  FSW-tool  pin 
diameter;  and  (d)  the  last  FSW  zone  which  is  located  above  the 
weld  nugget  is  generally  referred  to  as  the  flow  arm.  This  weld 
zone  typically  contains  the  workpiece  material  which  was 
(during  welding)  temporarily  confined  by  the  upright  truncated- 
conical  profile  of  the  FSW-tool  shoulder  underside. 

When  compared  to  the  conventional  fusion-welding  pro¬ 
cesses,  FSW  offers  a  number  of  advantages.  Since  a  fairly 
detailed  account  of  these  advantages  can  be  found  in  our  recent 
work  (Ref  11),  a  similar  detailed  account  of  these  advantages 
will  not  be  presented  here.  Nevertheless,  it  should  be  noted  that 
most  of  these  advantages  are  due  to  the  fact  that  FSW  is 
associated  with  lower  temperatures,  does  not  involve  fusion  and 
re-solidification  of  the  workpiece  material  and  that  no  filler 
metal,  flux  or  fuel/oxidizer  are  used.  Consequently,  the  extents 
of  microstructural  changes  and  the  associated  material  property 
degradation  are  smaller  relative  to  their  counterparts  observed 
in  the  conventional  fusion-based  welding  processes. 

Despite  the  fact  that  FSW  is  a  relatively  new  joining 
technology,  it  has  already  found  a  fairly  wide-scale  application 
in  several  industries.  Among  the  most  notable  examples  of  the 
successful  deployment  of  the  FSW  process  for  the  fabrication 
of  high-performance/cost-effective  structures  are:  (a)  Large- 
scale  production  of  aluminum  FSW-joined  ferryboat  deck 
structures  in  Finland;  (b)  manufacture  of  Al-Mg-Si-based  alloy 
FSW-joined  bullet-train  cabins  in  Japan;  (c)  fabrication  of 
Al-Cu-based  alloy  FSW-joined  rocket  launch  systems  by 
Boeing;  (d)  replacement  of  the  conventional  fusion  welding 
processes  with  FSW  in  the  manufacture  of  Al-Li  alloy-based 
space-shuttle  external  fuel-tanks  by  NASA;  and  (e)  application 
of  FSW  in  highly  demanding  jet  engine  applications  by  General 
Electric. 

Over  the  last  two  decades,  considerable  experimental 
research  efforts  have  been  invested  toward  obtaining  a  better 
understanding  of  the  FSW  joining  mechanism  and  the  accom¬ 
panying  evolution  of  the  welded-materials  microstructure/ 
properties  (e.g.,  Ref  12-15)  as  well  as  to  rationalizing  the 
effect  of  various  FSW  process  parameters  on  the  weld  quality/ 
integrity  (e.g.,  Ref  12,  16-19).  In  the  same  time  period, 
physically  based  modeling  and  simulations  of  the  FSW  process 
have  also  been  given  considerable  attention  (e.g.,  Ref  20-23).  A 


critical  assessment  of  the  past  experimental  and  computational 
approaches  revealed  that  the  majority  of  these  were  concerned 
with  important  practical  issues  such  as  the  development  of  tools 
and  identification  of  the  optimum  FSW-process  parameters  for 
a  variety  of  alloys  and  understanding  the  development  and 
characterization  of  the  weld-material  microstructure/properties. 
The  phenomenon  of  workpiece  material  flow,  despite  its  major 
role  in  the  fonnation,  structural  integrity,  and  the  overall 
properties  of  the  weld,  has  been  given  considerably  less 
attention.  As  a  main  objective  of  the  present  work,  an  attempt  is 
made  to  address  (computationally)  this  short  coming.  In 
particular,  it  is  argued  that  any  numerical  model  for  the  FSW- 
process  to  be  considered  reliable  and  of  high-fidelity  must 
correctly  predict  the  essential  features  of  the  material  flow  (as 
established  using  various  experimental  means,  section  2). 

The  organization  of  this  article  is  as  follows:  A  concise 
overview  of  the  major  experimental  FSW  flow-visualization 
techniques  and  findings  has  been  presented  in  section  2.  A  brief 
description  of  the  combined  Eulerian/Lagrangian  fully  coupled 
thermo-mechanical  model  and  computational  analysis  of  the 
FSW  process  is  provided  in  section  3.  The  key  computational 
results  pertaining  to  the  material  flow  during  the  FSW  process 
are  presented  and  compared  with  their  respective  experimental 
counterparts  in  section  4.  The  main  conclusions  resulting  from 
the  present  study  are  summarized  in  section  5. 


2.  Material  Flow  Visualization  During  FSW 

In  this  section,  a  brief  overview  is  provided  of  the  main 
FSW  flow  experimental-visualization  studies  carried  out  over 
the  past  decade.  One  of  the  first  FSW-flow  experimental- 
visualization  studies  was  conducted  by  Midling  (Ref  24)  who 
investigated  the  effect  of  the  FSW  tool  traverse  velocity  on  the 
workpiece  material  flow  (as  revealed  by  spatial  distribution  of 
the  material  interfaces  within  the  weld)  in  dissimilar-material 
FSW  joints.  This  procedure,  however,  provided  only  limited 
insight  into  the  workpiece  material  flow  since  it  relied  on  the 
visualization  of  the  dissimilar  alloys  interfaces  alone. 

Li  et  al.  (Ref  25)  carried  out  a  series  of  FSW  experiments 
involving  either  dissimilar  aluminum-alloy  grades  or  alumi¬ 
num/copper  workpiece  materials.  By  analyzing  material  pat¬ 
terns  seen  in  the  weld  metallographic  cross  sections  of  the  FSW 
welds,  Li  et  al.  established  that  the  resulting  workpiece  material 
flow  is  of  a  chaotic  nature. 

A  major  step  forward  in  the  visualization  of  the  material 
flow  was  made  in  the  work  of  Colligan  (Ref  26),  who  employed 
a  technique  in  which  small-diameter  steel  spheres  were  placed 
at  different  positions  along  the  workpiece  components  contact 
interface  before  welding  and  their  final  positions  within  the 
weld  were  determined  by  radiographing  the  post-welded 
structures.  The  main  findings  of  Colligan  (Ref  26)  can  be 
summarized  as  follows:  (a)  The  workpiece  material  is  stirred 
(via  chaotic-mixing)  only  in  the  upper  portion  of  the  weld;  and 
(b)  in  the  remaining  portions  of  the  weld  material  flow  involves 
only  extrusion  of  the  material  around  the  rotating  FSW  tool. 
The  main  limitations  of  Colligan’s  work  are:  (a)  FSW-flow 
visualization  was  inferred  by  tracking  only  single  points  (i.e., 
the  center  of  gravity  of  undeformed  steel  spheres);  and  (b)  the 
chaotic-mixing  observed  might  have  been  caused  by  the 
presence  of  finite-size  spheres  which  act  not  only  as  flow-field 
markers  but  also  cause  changes/disruptions  in  the  material  flow. 
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In  addition,  since  the  details  of  the  material  flow  are 
significantly  affected  by  the  FSW  tool  geometry  and  other 
FSW  process  parameters,  the  findings  of  Colligan  (Ref  26) 
could  not  be  readily  generalized. 

The  next  major  step  forward  in  the  visualization  of  the 
material  flow  was  made  in  a  series  of  papers  by  Reynolds  and 
co-workers  (Ref  27-30)  who  introduced  the  so-called  marker 
insert  technique.  Within  this  technique,  markers  of  a  different 
aluminum-alloy  grade  are  inserted  along  and  on  both  sides  of 
an  aluminum-alloy  workpiece  contacting  interface  (i.e.,  in  the 
path  of  a  rotating  and  traversing  FSW  tool).  The  final  position 
and  shape  of  the  markers  were  determined  by:  (a)  Milling-off 
successive  0.25  mm  thick  slices  from  the  top  surface  of  the 
weldment;  (b)  etching  each  newly  created  top  surface  of  the 
weldment  using  an  etchant  which  produces  a  strong  contrast 
between  the  two  aluminum-alloy  grades;  (c)  capturing  a  digital 
image  of  each  etched  weldment  top  surface;  (d)  digital  image 
analysis  of  each  etched  weldment  top  surface;  and  (e)  assembly 
of  the  individual  images  to  fonn  a  complete  three-dimensional 
model  for  the  FSW  weldment  revealing  the  position  and  the 
deformed  shape  of  the  markers.  These  findings  are  subse¬ 
quently  used  to  infer  the  essential  features  of  material  flow 
during  the  FSW  process. 

The  main  findings  obtained  by  Reynolds  and  co-workers 
(Ref  27-30)  can  be  summarized  as  follows:  (a)  The  material 
flow  during  FSW  is  influenced  by  both  the  translational  and 
rotational  components  of  the  FSW  tool  motion;  (b)  the  main 
component  of  the  material  flow  is  in  the  horizontal  plane  and  in 
a  general  direction  which  is  opposite  to  the  direction  of  FSW 
tool  travel  and  involves  material  motion/extrusion  around  the 
tool  and  its  placement  in  a  region  behind  the  tool.  In  this  case, 
the  tool  shoulder,  the  pin,  the  backing  plate,  and  the  base-metal 
portion  of  the  workpiece  material  effectively  form  an  extrusion 
die;  (c)  workpiece  material  stirring  (i.e.,  its  rotational  motion 
around  the  tool  longitudinal  axis)  occurs  predominantly  in  the 
top  portion  of  the  weld  in  which  the  rotating  action  of  the  tool 
shoulder  causes  material  transport  from  the  retreating  to  the 
advancing  side  of  the  weld;  (d)  the  material  transferred  from  the 
retreating  to  the  advancing  side  interferes  with  the  extrusion  (in 
the  horizontal  plane)  of  the  advancing-side  material  in  the  path 
of  FSW-tool  pin.  This  results  in  a  downward  (out-of-plane) 
component  of  the  material  extrusion  direction.  The  downward 
motion  of  the  extruded  material  extends  through  the  thickness 
of  the  workpiece  material  on  the  advancing  side  till  the  weld 
root  at  which  point,  due  to  the  presence  of  the  rigid  backing 
plate,  the  extruded  material  crosses  over  to  the  retreating  side 
and  begins  its  ascent  (i.e.,  to  move  in  the  upward  direction).  As 
a  consequence,  the  material  acquires  a  rotational  motion  around 
the  longitudinal  axis  of  the  weld;  (e)  the  extent  of  the  vertical 
circular  motion  described  in  (d)  increases  with  a  reduction  in 
the  value  of  the  weld-pitch  which  is  attributed  to  the  auguring 
effect  of  the  threaded  FSW-tool  pin;  and  (f)  the  extent  of 
vertical  rotation  also  increases  with  the  pin  diameter  and  this 
effect  is  attributed  to  the  associated  increase  in  the  amount  of 
energy  transferred  into  the  weld  by  the  FSW  tool. 


3.  FSW  Computational  Modeling  and  Analysis 

Modeling  of  the  FSW  process  was  carried  out  in  our  prior 
work  using  an  Arbitrary  Lagrangian  Eulerian  (ALE)  fully 
coupled  thermo-mechanical  finite-element  procedure  (Ref  20- 


22).  The  main  emphasis  of  our  prior  work  was  the  prediction  of 
microstructure  and  property  fields  within  the  weld  region  in  a 
variety  of  solid-solution  strengthened  and  age-hardened  alumi¬ 
num-  and  titanium-based  alloys.  Since  the  main  emphasis  of  the 
present  work  is  the  analysis  of  the  material  flow  during  the  FSW 
process  under  different  tool-design  and  process-parameter 
conditions,  a  new  Combined  Eulerian  Lagrangian  (CEL)  fully 
coupled  thermo-mechanical  analysis  of  the  FSW  process  is 
introduced.  While  this  analysis  is  usually  computationally  more 
expensive  than  the  corresponding  ALE  analysis,  it  is  generally 
found  to  be  more  robust  and  more  suitable  for  studying  material¬ 
mixing  processes.  It  should  be  recognized,  however,  that  the 
ALE  and  CEL  analyses  share  many  features.  Flence,  the  features 
common  to  the  two  analyses  will  only  be  discussed  briefly  here. 

3.1  FSW  Computational  Analysis 

3.1.1  Computational  Domain.  The  computational 
domain  used  consists  of  two  separate  sub-domains,  one  of  an 
Eulerian-type  and  the  other  of  a  Lagrangian-type.  The  Eulerian 
sub-domain  (used  to  model  the  workpiece)  is  of  a  parallelepiped 
shape  with  the  dimensions  of  L  x  Wx  H  =  50x  40  x8  mm, 
respectively.  This  sub-domain  is  typically  discretized  using 
Eulerian  8-node  brick  elements  with  a  characteristic  edge  length 
of  0.44  mm.  The  Lagrangian  sub-domain  is  used  to  model  the 
FSW  tool.  As  will  be  discussed  in  section  4,  a  prototypical  FSW 
tool  design  consisting  of  a  tapered  threaded  pin  and  a  cylindrical 
shoulder  with  an  upright  truncated  conical  under-cut  was 
utilized  in  the  present  work.  The  tools  are  meshed  using 
4-node  tetrahedron  continuum  elements.  Examples  of  the 
Eulerian  and  Lagrangian  sub-domains  used  in  the  present  work 
are  displayed  in  Fig.  4.  It  should  be  noted  that  in  this  figure,  for 
clarity,  the  tool  is  shown  in  the  retracted  position.  In  the  present 
work,  the  Eulerian  sub-domain  typically  consisted  of  150,000 
elements  while  the  tool  contained  20,000  elements. 

3.1.2  Computational  Analysis  Type.  The  FSW  process 
is  analyzed  computationally  using  a  Combined  Eulerian 
Lagrangian  and  a  fully  coupled  thermo-mechanical  finite- 
element  algorithm.  Within  this  algorithm,  heat  dissipation 
associated  with  plastic  defonnation  and  tool/workpiece  inter¬ 
facial  friction-sliding  is  treated  as  a  heat-source  in  the 
governing  thermal  equation.  On  the  other  hand,  the  effect  of 
temperature  on  the  mechanical  response  of  the  workpiece 
material  is  taken  into  account  through  the  use  of  a  temperature- 
dependent  workpiece  material  model. 

3.1.3  Initial  Conditions.  The  Eulerian  sub-domain  is 
initially  filled  with  the  workpiece  and  the  flow-visualization 
marker  materials  by  prescribing  the  appropriate  material 


Fig.  4  A  schematic  of  a  typical  FSW  tool  (Lagrangian  domain)  and 
workpiece  (Eulerian  domain)  used  in  the  present  work 
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volume  fractions  to  each  Eulerian  element.  To  enable  the 
motions  of  the  Eulerian  materials  at  the  workpiece  upper 
surface  in  the  outward  normal  direction  (without  material  loss), 
one  or  more  top  Eulerian  element  layers  is  initially  left  void. 
For  the  Lagrangian  sub-domain,  at  the  beginning  of  the  analysis 
the  tool  is  assigned  a  constant  rotational  velocity  and  a  zero 
translational  velocity.  Both  the  Eulerian  and  Lagrangian  sub- 
domains  are  assigned  ambient-temperature  initial  conditions. 

3.1.4  Boundary  Conditions.  For  convenience,  the  longi¬ 
tudinal  motion  of  the  FSW  tool  is  not  considered  explicitly. 
Instead,  the  longitudinal  velocity  of  the  tool  is  set  to  zero  and  the 
workpiece  material  is  allowed  to  move  through  the  Eulerian  sub- 
domain  with  an  overall  longitudinal  velocity  equal  to  negative 
of  the  tool  translational  velocity.  This  was  accomplished  through 
the  use  of  the  appropriate  “in-flow”  and  “out-flow”  velocity 
boundary  conditions  over  the  vertical  faces  of  the  Eulerian 
domains  which  are  orthogonal  to  the  direction  of  tool  travel. 
Thus,  the  Eulerian  sub-domain  displayed  in  Fig.  4  does  not 
represent  the  entire  workpiece  but  rather  a  rectangular  region 
around  the  tool  in  the  otherwise  infinitely-long  workpiece. 

“No-flow”  boundary  conditions  are  prescribed  along  the 
Eulerian  sub-domain  vertical  faces  parallel  to  the  tool  travel 
direction  to  account  for  the  role  of  the  workpiece  base-metal 
material  in  restraining  material  flow  in  the  transverse  direction. 
To  mimic  the  role  played  by  the  workpiece  rigid  backing  plate 
in  preventing  the  flow  of  the  workpiece  material  in  the 
downward  direction,  zero  normal  velocity  boundary  conditions 
are  applied  over  the  bottom  surface  of  the  Eulerian  sub-domain. 
In  accordance  with  the  initial  conditions  prescribed  to  the  top 
Eulerian-element  layers,  “out-flow”  boundary  conditions  are 
prescribed  over  the  top  surface  of  the  Eulerian  sub-domain. 

As  far  as  the  thermal  boundary  conditions  are  concerned, 
standard  convective  boundary  conditions  are  applied  over  free 
surfaces  of  the  workpiece  and  the  tool  while  enhanced 
convection  boundary  conditions  are  applied  over  the  bottom 
face  of  the  workpiece  (to  mimic  the  effect  of  enhanced  heat 
extraction  through  the  workpiece  backing  plate). 

3.1.5  Tool/Workpiece  Contact  Interactions.  Since  the 
Eulerian  and  Lagrangian  domains  do  not  possess  conformal 
meshes,  the  contact  interfaces  between  the  two  could  not  be 
defined  using  mesh-based  surfaces.  Instead,  contact  interfaces 
between  the  Lagrangian  and  the  Eulerian  sub-domains  are 
determined  using  the  so-called  “immersed  boundary  method” 


(Ref  31)  which  identifies,  during  each  computational  time 
increment,  the  boundary  of  the  Eulerian  sub-domain  region 
which  is  occupied  by  the  Lagrangian  sub-domain.  Eulerian- 
Lagrangian  contact  constraints  are  enforced  using  a  penalty 
method,  within  which  the  extent  of  contact  pressure  is  governed 
by  the  local  surface  penetrations  (where  the  default  penalty 
stiffness  parameter  is  automatically  maximized  subject  to 
stability  limits).  As  far  as  the  shear  stresses  are  concerned  they 
are  transferred  via  a  “slip/stick”  algorithm  that  is  shear  stresses 
lower  than  the  frictional  shear  stress  are  transferred  without 
interface  sliding  (otherwise  interface  sliding  takes  place).  The 
frictional  shear  stress  is  defined  by  a  modified  Coulomb  law 
within  which  there  is  an  upper  limit  to  this  quantity  (set  equal  to 
the  shear  strength  of  the  workpiece  material).  The  frictional  shear 
stress  is  then  defined  as  a  smaller  of  the  product  between  the 
static/kinetic  friction  coefficient  and  the  contact  pressure,  on  one 
hand,  and  the  workpiece  material  shear  strength,  on  the  other. 

In  addition  to  the  Eulerian-Lagrangian  contacts,  interactions 
(of  a  “sticky”  character)  also  occur  between  different  Eulerian 
materials.  This  type  of  interaction  is  a  consequence  of  the 
kinematic  constraint  which  requires  that  all  Eulerian  materials 
residing  in  a  single  Eulerian  element  are  subjected  to  the  same 
strain.  The  Eulerian-Eulerian  contacts  allow  nonnal  (tensile  and 
compressive)  stresses  to  be  transferred  between  adjoining 
materials  while  no  slip  at  the  associated  material  boundaries  is 
allowed. 

3.1.6  Heat-Generation  and  Partitioning.  As  mentioned 
earlier,  both  plastic  defonnation  and  frictional  sliding  are  treated 
as  heat  sources.  To  account  for  the  fact  that  a  small  fraction  of 
the  plastic -deformation  work  is  stored  in  the  form  of  crystal 
defects,  95%  of  this  work  was  assumed  to  be  dissipated  in  the 
form  of  heat.  As  far  as  heat  generation  due  to  frictional  sliding  is 
concerned,  it  is  assumed  that  its  rate  scales  with  the  product  of 
local  interfacial  shear  stress  and  the  sliding  rate,  and  that  100% 
of  this  energy  is  dissipated  in  the  form  of  heat.  Partitioning  of 
this  heat  between  the  tool  and  the  workpiece  is  then  computed 
using  the  algorithm  reported  in  Ref  32  and  the  appropriate 
thermal  properties  of  the  workpiece  and  tool  materials. 

3.1.7  Computational  Algorithm.  As  mentioned  earlier,  a 
new  CEL-based  finite  element  analysis  of  the  FSW  process  is 
developed.  Within  this  analysis,  the  workpiece  is  treated  as  an 
Eulerian  sub-domain,  the  tool  is  treated  as  a  Lagrangian  sub- 
domain  and  the  interaction  between  the  two  is  treated  using  an 


Table  1  Johnson-Cook  yield  strength  (<Ty)  material  model  parameters*  and  the  corresponding  elastic  and  thermal 
parameters  for  AA5059-H131 


Parameter 

Symbol 

Units 

Value 

Reference  strength 

A 

MPa 

167.0 

Strain-hardening  parameter 

B 

MPa 

596.0 

Strain-hardening  exponent 

n 

N/A 

0.551 

Strain-rate  coefficient 

c 

N/A 

0.001 

Room  temperature 

T 

±  room 

K 

293 

Melting  temperature 

Tmelt 

K 

893.0 

Temperature  exponent 

in 

N/A 

1.0 

Young’s  modulus 

E 

GPa 

70 

Poisson’s  ratio 

V 

N/A 

0.3 

Density 

p 

kg/m3 

2700 

Thermal  conductivity 

k 

W/m-  K 

120 

Specific  heat 

cp 

J/kg-K 

880 

=  [a  -hS^1)"]  [1  +  ciog^'/lc1)]  [1  -  Tff];  TH  =  {t  -  rroom)/(rmclt  -  rroom) 
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Fig.  5  Side  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  advancing-side/ 
top-layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e) 
18.3,  (f)  22.2,  (g)  26.1,  (h)  30 


Eulerian/Lagrangian  contact  (“fluid-structure  interaction”) 
algorithm.  Within  the  Lagrangian  sub-domain:  (a)  The  mesh 
(nodes  and  elements)  is  attached  to  the  associated  material  and 
moves  and  deforms  with  it;  and  (b)  each  element  must  be  fully 
filled  with  a  single  material.  On  the  other  hand,  within  an 
Eulerian  sub-domain:  (a)  The  mesh  is  fixed  in  space  and  the 
material  flows  through  it;  (b)  elements  are  allowed  to  be  partially 
filled  and/or  contain  multiple  materials;  and  (c)  since  the 
material  and  the  element  boundaries  do  not  generally  coincide,  a 
separate  (“interface  reconstruction”)  algorithm  must  be  used  to 
track  the  position  of  Eulerian  material  boundaries.  The  interface 
reconstruction  algorithm  approximates  the  material  boundaries 


within  an  element  as  simple  planar  facets  and,  hence,  accurate 
determination  of  a  material’s  location  within  an  element  requires 
the  use  of  fine  Eulerian  meshes. 

Numerical  solution  of  the  governing  equations  in  the  Eulerian 
sub-domain  within  each  time  increment  involves  two  separate 
steps:  (a)  the  Lagrangian  step  within  which  the  sub-domain  is 
temporarily  treated  as  being  of  a  Lagrangian-type  (i.e.,  nodes  and 
elements  are  attached  to  and  move/deform  with  the  material); 
and  (b)  the  “remap”  step  within  which  the  distorted  mesh  is 
mapped  onto  the  original  Eulerian  mesh  and  the  accompanying 
material  transport  is  computed  and  used  to  update  the  Eulerian- 
material  states  and  inter-material  boundaries. 
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Fig.  6  Top  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  advancing-side/top- 
layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e) 
18.3,  (f)  22.2,  (g)  26.1,  (h)  30 


It  should  be  noted  that  the  main  reasons  for  using  the  CEL 
approach  is  that  the  conventional  Lagrangian  finite  element 
analysis:  (a)  Does  not  allow  material  mixing  within  individual 
finite  elements;  and  (b)  becomes  inaccurate  and  failure-prone 
when  the  elements  experience  a  large  degree  of  distortion. 

The  fully  coupled  thermo-mechanical  problem  dealing  with 
FSW  is  solved  using  an  explicit  solution  algorithm  implemented  in 
ABAQUS/Explicit  (Ref  3 1),  a  general  purpose  finite  element  solver. 

3.1.8  Computational  Accuracy,  Stability,  and  Cost.  A 
standard  mesh  sensitivity  analysis  was  carried  out  (the  results 
not  shown  for  brevity)  to  ensure  that  the  results  obtained  are 
accurate,  i.e.,  insensitive  to  the  size  of  the  elements  used. 


Due  to  the  conditionally  stable  nature  of  the  explicit  finite 
element  analysis  used,  the  maximum  time  increment  during 
each  computational  step  had  to  be  lower  than  the  attendant 
stable  time  increment.  To  keep  the  computational  cost 
reasonable  while  ensuring  accuracy  and  stability  of  the 
computational  procedure,  a  mass  scaling  algorithm  is  used. 
This  algorithm  adaptively  adjusts  material  density  in  the 
critical  stable  time-increment  controlling  finite  elements  with¬ 
out  significantly  affecting  the  computational  analysis  results.  A 
typical  30  s  FSW  computational  analysis  required  10-h  of 
(wall-clock)  time  on  a  12  core,  3.0  GHz  machine  with  16  GB 
of  memory. 
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Fig.  7  Rear  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  advancing-side/ 
top-layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e) 
18.3,  (f)  22.2,  (g)  26.1,  (h)  30 


3.2  Material  Models 

3.2.1  Tool  Material.  FSW  tools  used  for  joining  alumi¬ 
num  alloys  are  typically  made  of  one  of  the  hot-working  or  high¬ 
speed  steel  grades  (e.g.,  AISI  H13  or  AISI  350M).  Due  to  high 
strength  of  these  steel  grades  at  elevated  temperatures,  the  tool 
typically  experiences  very  little  plastic  deformation  during  the 
FSW  process.  On  the  other  hand,  tool  wear  and  loss  of  tool 
features  after  prolonged  use  is  frequently  observed.  Since  an 
analysis  of  tool  wear  is  beyond  the  scope  of  the  present  work 
and  the  likely  hood  for  plastic  deformation  of  the  tool  is  very 
small,  the  tool  material  (AISI  HI  3)  is  modeled  as  a  mechanically 


isotropic  linear-elastic  material  with  a  Young’s  Modulus  of 
210  GPa,  a  Poisson’s  ratio  of  0.3,  and  density  of  7825  kg/m3. 
The  thermal  properties  of  this  material  which  play  an  important 
role  in  partitioning  of  frictional-sliding-induced  heat  at  the  tool/ 
workpiece  interface  are  assigned  as:  thermal  conductivity, 
k  =  28.5  W/m  •  K  and  specific  heat,  cp  =  475  J/kg  •  K. 

3.2.2  Workpiece  Material.  The  workpiece  material 
(AA5059-H131,  a  solid-solution  strengthened  and  strain- 
hardened/stabilized  Al-Mg-Mn  alloy)  is  assumed  to  be  isotro¬ 
pic,  linear-elastic  and  strain-hardenable,  strain-rate  sensitive, 
thermally  softenable  plastic  material  and  is  modeled  using  the 
Johnson-Cook  material  model  (Ref  33).  A  summary  of  the 
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Fig.  8  Side  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  interface/top-layer) 
volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt  angle  =  2.5°; 
tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e)  18.3,  (f)  22.2,  (g) 
26.1,  (h)  30 


Johnson-Cook  material  model  parameters  (including  the  yield 
stress,  <Jy,  governing  equation)  and  the  elastic  and  thermal 
properties  of  AA5059-H131  are  provided  in  Table  1. 

Examination  of  the  original  Johnson-Cook  material  model 
reveals  that  a  temperature  change  causes  only  a  reversible 
change  in  material  strength  by  affecting  thennal  activation  of 
dislocation  motion.  In  other  words,  the  original  Johnson-Cook 
model  does  not  account  for  any  pennanent  changes  in  the 
material  microstructure  and  properties  which  may  occur  as  a 
result  of  high-temperature  exposure  of  the  material.  This 
approach  is  not  fully  justified  in  the  case  of  FSW  where,  due  to 
the  attendant  high  temperatures,  significant  differences  in  the 
material  microstructure  and  properties  may  exist  between  the 


base  metal  and  the  weld  (as  well  as  within  different  regions  of 
the  weld).  To  account  for  this  additional  effect  of  temperature,  a 
modification  of  the  strain  hardening  term  within  the  original 
Johnson-Cook  model  was  proposed  in  our  prior  work  (Ref  20). 
Specifically,  the  strain  hardening  term  is  still  assumed  to  be  a 
parabolic  function  of  equivalent  plastic  strain  (  Be!pX,  where  B 
and  n  are  material  parameters,  Table  1).  However,  the  equiv¬ 
alent  plastic  strain,  epi,  is  now  defined  as  the  sum  of  two  terms: 
one  (positive)  which  quantifies  the  contribution  of  plastic 
deformation  to  strain  hardening  and  the  other  (negative)  which 
accounts  for  strain  softening  induced  by  dynamic  recrystalli¬ 
zation.  In  other  words,  strain  hardening  is  still  assumed  to  be 
controlled  by  the  density  of  (mobile)  dislocations  but  the  local 
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Fig.  9  Top  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  interface/top-layer) 
volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt  angle  =  2.5°; 
tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e)  18.3,  (f)  22.2,  (g) 
26.1,  (h)  30 


544  value  of  this  quantity  is  taken  to  be  the  result  of  a  competition 

545  between  dislocation-generating  plastic  deformation  and  dislo- 

546  cation  annihilation  associated  with  dynamic  recrystallization.  It 

547  should  be  noted  that  the  first  (positive)  component  of  the 

548  equivalent  plastic  strain  in  the  modified  Johnson-Cook  material 

549  model  still  quantifies  the  overall  extent  of  inelastic  deformation 

550  taking  place  at  a  given  material  point. 

551  4.  Results  and  Discussion 

552  As  mentioned  earlier,  the  main  objective  of  the  present  work 

553  is  to  study  the  flow  of  the  workpiece  material  during  FSW  and 


to  correlate  this  flow  with  the  values/levels  of  the  friction  stir 
welding  process  key  parameters.  Toward  that  end,  a  marker  was 
placed  (on  either  the  advancing  or  retreating  workpiece  sides, 
or  along  the  weld  interface)  along  the  path  of  the  advancing 
FSW  tool.  The  marker  material  was  identical  to  the  remainder 
of  the  workpiece  material.  Flowever,  by  declaring  the  marker 
material  as  a  different  material  within  the  Eulerian  region,  it 
was  possible  to  monitor  the  flow  of  the  workpiece  material 
during  the  FSW  process  and  to  assess  the  extent  of  material 
mixing  during  FSW. 

To  help  with  the  visualization  of  the  marker  material  flow,  the 
volume  fraction  of  the  marker  material  within  the  Eulerian  cells 
was  monitored.  At  the  beginning  of  the  simulation,  each  marker 
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Fig.  10  Rear  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  interface/top- 
layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  2.7,  (b)  6.6,  (c)  10.5,  (d)  14.4,  (e) 
18.3,  (f)  22.2,  (g)  26.1,  (h)  30 


(one  marker  was  used  in  a  single  analysis)  contained  hundred 
percent  of  the  marker  material  (i.e.,  the  volume  fraction  of  the 
marker  material  within  the  marker  was  initially  1.0).  However, 
as  the  friction  stir  welding  process  simulation  proceeded,  due  to 
the  flow  of  the  marker  material  and  its  mixing  with  the 
remainder  of  the  workpiece  material,  the  marker  material  spread 
over  a  wider  region.  At  the  same  time,  in  accordance  with  the 
mass  conservation  requirement,  as  the  marker  material  was 
spread  over  a  wider  region  the  volume  fraction  of  the  marker 
material  experienced  a  continues  decrease. 

To  further  help  visualize  the  material  flow  during  FSW,  the 
workpiece  material  is  not  shown  in  the  figures  displayed  in  this 


section.  Only  the  marker  material  distribution/volume  fraction 
is  displayed  (in  addition  to  that  of  the  friction  stir  welding  tool). 
Also,  the  friction  stir  welding  tool  was  made  partially 
transparent  so  that  the  location  of  the  marker  material  under 
the  friction  stir  welding  tool  shoulder  can  be  monitored. 

In  this  section,  the  results  of  the  material  flow  analysis  are 
presented  and  discussed.  The  results  are  presented  in  the 
following  way:  First,  a  base  line  case  is  considered  within 
which  a  prototypical  set  of  friction  stir  welding  process 
parameters  was  chosen.  This  is  followed  by  a  section  in  which 
the  effect  of  variation  of  few  specific  FSW  process  parameters 
was  investigated. 


Journal  of  Materials  Engineering  and  Perfonnance 


Journal  :  JMEP 

11665 

Dispatch:  24-10-2011 

□  LE 

Pages  :  17 

□  TYPESET 

Article  Number : 

69 

lif  CP 

0  DISK 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 

601 


602 

603 

604 

605 

606 

607 

608 

609 

610 

611 

612 


Fig.  11  Side  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  retreating  side/ 
top-layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  3.9,  (b)  5.4,  (c)  7.2,  (d)  8.4,  (e)  9.9, 
(f)  11.4,  (g)  12.9,  (h)  14.4 


4.1  Baseline  Case 

In  this  section,  prototypical  values/levels  are  assigned  to  the 
FSW  key  process  parameters.  Specifically:  (a)  tool  geometry:  (i) 
tool  shoulder  diameter  =  0.025  m,  (ii)  FSW  tool  shoulder 
truncated  conical  under-cut  angle  =  7°,  (iii)  pin  length  = 
0.0057  m,  (iv)  pin  upper  diameter  =  0.01  m,  (v)  pin-taper 
angle  =  10°,  and  (vi)  pin  thread  pitch  =  800  threads/m;  (b)  tool 
material  =  AISI  FI  13  tool  steel;  (c)  workpiece  material  = 
AA5059;  (d)  tool  rotation  speed  =  500  rpm;  (e)  tool  travel 
speed  =  0.0025  m/s;  (f)  tool  tilt  angle  =  2.5°,  and  (g)  tool  plunge 
depth  =  0.0001  m. 


Spatial  distribution  and  temporal  evolution  of  the  volume 
fraction  of  the  marker  material  for  the  base  line  case  in  which  the 
marker  was  placed  on  the  advancing  side  of  the  workpiece  are 
depicted  in  Fig.  5(a)  to  (h).  The  results  displayed  in  these  figures 
represent  a  side  view  in  which  the  tool  is  traveling  to  the  left 
(relative  to  the  workpiece  material).  The  corresponding  top  and 
rear  (behind  the  tool)  views  of  the  same  results  are,  respectively, 
displayed  in  Fig.  6(a)  to  (h)  and  7(a)  to  (h),  respectively. 

The  corresponding  baseline  case  results  in  the  three  views  for 
the  case  when  the  marker  was  placed  at  the  welding  interface  are 
displayed  in  Fig.  8(a)  to  (h),  9(a)  to  (h),  and  10(a)  to  (h). 
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Fig.  12  Top  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  retreating  side/ 
top-layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  3.9,  (b)  5.4,  (c)  7.2,  (d)  8.4,  (e)  9.9, 
(f)  11.4,  (g)  12.9,  (h)  14.4 


Likewise,  the  corresponding  results  for  the  case  when  the  marker 
was  placed  on  the  retreating  side  of  the  workpiece  are  shown  in 
Fig.  11(a)  to  (h),  12(a)  to  (h),  and  13(a)  to  (h),  respectively. 

It  should  be  noted  that  the  results  presented  in  Fig.  5(a)  to 
(h)  to  13(a)  to  (h)  correspond  to  the  cases  when  the  marker  was 
placed  at  the  top  surface  of  the  workpiece  (as  labeled  in  part  (a) 
of  these  figures).  Within  the  present  work,  the  workpiece 
configurations  in  which  the  marker  was  placed  half  way 
through  the  workpiece  thickness  or  near  the  bottom-surface  of 
the  workpiece  were  also  analyzed.  Since  these  results  were 
qualitatively  similar  to  the  ones  presented  here  (although  they 
showed  a  decrease  in  the  extent  of  material  mixing),  they  will 
not  be  presented  here. 


Examination  of  the  results  presented  in  Fig.  5(a)  to  (h)  to 

13(a)  to  (h)  revealed  that: 

(a)  While  the  marker  material  was  initially  confined  within  a 
cubic  domain  (the  initial  configuration  is  not  shown  for 
brevity),  due  to  its  flow  and  mixing  with  the  surrounding 
workpiece  material,  it  became  dispersed  over  a  larger 
spatial  domain  of  the  workpiece.  At  the  same  time,  the 
volume  fraction  of  the  marker  material  decreased,  in 
accordance  with  the  mass  conservation  requirements; 

(b)  At  the  completion  of  the  FSW  process,  the  marker  mate¬ 
rial  tends  to  reside  on  the  same  side  of  the  welding 
interface  as  that  where  it  was  originally  placed  (e.g.,  the 
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Fig.  13  Rear  view  of  the  spatial  distribution  and  temporal  evolution  of  the  marker-material  (initially  located  on  the  workpiece  retreating  side/ 
top-layer)  volume  fraction  (red  greater  than  0.8,  blue  less  than  0.2):  tool  rotation  speed  =  500  rpm;  tool  travel  speed  =  0.0025  m/s;  tool  tilt 
angle  =  2.5°;  tool  plunge  depth  =  0.0001  m;  10  mm  upper-diameter  tool  pin;  welding  time  in  seconds:  (a)  3.9,  (b)  5.4,  (c)  7.2,  (d)  8.4,  (e)  9.9, 
(f)  11.4,  (g)  12.9,  (h)  14.4 


advancing  side  marker-material  is  ultimately  placed  in 
the  advancing  side  weld  region).  However,  it  should  be 
noted  that  this  process  is  accompanied  by  a  considerable 
dispersion  of  the  marker  material  in  both  the  lateral  and 
through-the-thickness  directions; 

(c)  For  the  most  part,  FSW  causes  the  workpiece  material  flow 
in  the  horizontal  plane  (a  plane  parallel  with  the  workpiece 
backing  plate),  e.g.  Fig.  5(b)  to  (d).  In  addition  (down¬ 
ward)  flow  of  the  marker  material  in  the  through-the-thick- 
ness  direction  was  also  observed.  In  the  case  of  the 
advancing-side  marker,  this  component  of  the  material 
flow  was  fairly  pronounced,  Fig.  5(h)  and  7(h).  This  com¬ 


ponent  of  the  material  flow  still  exists  but  is  less  pro¬ 
nounced  in  the  case  of  the  weld-interface  marker,  Fig.  8(h) 
and  10(h).  On  the  other  hand,  practically  no  marker  mate¬ 
rial  flow  in  the  through-the-thickness  direction  occurs  in 
the  case  of  the  retreating-side  marker,  Fig.  11(h)  and  13(h); 

(d)  The  material  downward  flow  is,  at  least  partly,  promoted 
by  the  threaded-pin  auguring  effect,  e.g.,  Fig.  8(f)  and  (g); 

(e)  The  extent  of  marker-material  dispersion  is  highest  in  the 
case  of  the  advancing  side  marker,  Fig.  6(h).  It  is  inter¬ 
mediate  in  the  case  of  the  weld-interface  marker  and  the 
lowest  in  the  case  of  the  retreating-side  marker.  In  other 
words,  the  advancing  material  undergoes  extensive 
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Fig.  14  Rear  view  of  the  spatial  distribution  of  the  marker-material  (initially  located  on  the  workpiece  advancing-side/top-layer)  volume  frac¬ 
tion  for:  (a)  baseline  case,  tool  rotation  speed  =  500  rpm,  tool  travel  speed  =  0.0025  m/s,  tool  plunge  depth  =  0.0001  m,  tool  tilt  angle  =  2.5°, 
10  mm  upper-diameter  tool  pin;  (b)  tool  rotation  speed  =  700  rpm;  (c)  tool  rotation  speed  =  300  rpm;  (d)  tool  tilt  angle  =  3.50°;  (e)  tool  tilt 
angle  =  2.0°;  (f)  12  mm  upper-diameter  tool  pin;  and  (g)  8  mm  upper-diameter  tool  pin.  Note  that  for  cases  (b)  to  (g)  the  values  of  the  unspeci¬ 
fied  FSW  parameters  are  the  same  as  in  the  baseline  case  (a) 


stirring  before  it  is  ultimately  placed  into  the  weld.  In 
sharp  contrast,  the  retreating-side  marker  material  is 
mainly  extruded  around  the  tool  and  forged  in  the  wake  of 
the  tool  (without  experiencing  substantial  mixing);  and 
(f)  There  is  no  evidence  for  rotation  of  the  marker  material 
around  the  welding  direction  in  any  of  the  results  being 
analyzed.  That  is,  no  instances  of  marker  material  cross¬ 
ing  from  the  advancing  to  the  retreating  side  were 
observed  at  the  lower  portion  of  the  workpiece  and  in 
the  wake  of  the  FSW  tool  as  suggested  by  the  experi¬ 
mental  results  of  Reynolds  and  co-workers  (Ref  27-30). 


4.2  The  Effect  of  FSW-Process  Parameters  on  Workpiece 
Material  Flow 

Within  the  present  work,  the  effect  of  the  following  FSW 
process  parameters  on  the  workpiece  material  flow  has  been 
investigated:  (a)  Weld  pitch  (defined  as  the  ratio  of  the  FSW  tool 
travel  speed  to  the  tool  rotation  speed);  (b)  tool  tilt-angle;  and  (c) 
the  tool  pin  upper  and  lower  diameters  (at  a  constant  pin  taper 
angle).  In  all  the  simulations  presented  in  this  section,  only  one 
of  these  process  parameters  was  altered  (i.e.,  all  the  remaining 
process  parameters  were  kept  unchanged  relative  to  their 
baseline  counterparts).  To  help  in  revealing  the  effect  of  these 
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process  parameters  on  the  workpiece  material  flow,  the  rear¬ 
view,  end-of-the-simulation,  advancing-side,  marker-material 
volume  fraction  results  for  the  base-line  case  are  displayed  in 
Fig.  14(a).  A  comparison  of  the  results  displayed  in  Fig.  14(a) 
with  the  results  displayed  in  the  remaining  parts  of  Fig.  14  are 
used  to  infer  the  effect  of  the  specific  FSW  process  parameters. 

4.2.1  Weld  Pitch.  The  effect  of  the  weld  pitch  can  be 
inferred  by  comparing  the  results  displayed  in  Fig.  14(b),  (a), 
and  (c),  respectively.  The  magnitudes  of  the  weld  pitch  in  these 
three  cases  are:  3.41  x  10~5  m/radian,  4.77  x  1(P5  m/radian, 
and  7.95  x  10~5  m/radian,  respectively.  Examination  of  the 
results  displayed  in  these  figures  reveals  that  the  highest  extent 
of  marker-material  mixing/dispersion  is  attained  at  the  lowest 
value  of  the  weld  pitch.  Fig.  14(b)  and  this  extent  continuously 
decreases  with  an  increase  in  the  magnitude  of  the  weld  pitch. 
Fig.  14(a)  and  (c).  This  finding  is  expected  since  a  lower  value 
of  the  weld  pitch  (i.e.,  an  increased  tool  rotational  speed  at  a 
constant  tool  travel  speed)  represents  an  increase  in  the 
mechanical  energy  introduced  into  the  workpiece  through  the 
rotating  FSW  tool.  Dissipation  of  this  energy  in  the  form  of 
heat  increases  the  workpiece  temperature,  softens  the  work- 
piece  material  and,  hence,  facilitates  material  mixing. 

4.2.2  Tool  Tilt-Angle.  The  effect  of  the  tool  tilt-angle  can 
be  inferred  by  comparing  the  results  displayed  in  Fig.  14(d), 
(a),  and  (e),  respectively.  The  magnitudes  of  the  tool  tilt-angle 
in  these  three  cases  are:  3.5°,  2.5°,  and  2.0°,  respectively. 
Examination  of  the  results  displayed  in  these  figures  reveals 
that  the  highest  extent  of  marker-material  mixing/dispersion  is 
attained  in  the  baseline  case,  Fig.  14(a)  and  that  this  extent 
decreases  with  a  deviation  in  the  tilt-angle  from  its  baseline 
value.  In  other  words,  a  tilt-angle  around  the  baseline  value  of 
2.5°  appears  to  be  optimal.  This  finding  is  in  accordance  with 
the  experimental  observations  in  the  FSW  welding  practice. 
The  existence  of  an  optimum  tilt-angle  is  generally  described  as 
a  trade-off  between  the  maximization  of  the  forging  pressure 
(favors  a  large  tilt-angle)  and  minimization  of  the  workpiece- 
weld  “ploughing”  (favors  a  smaller  tilt-angle).  The  present 
results  show  that  maximization  of  workpiece-material  stirring/ 
mixing  also  favors  an  optimal  value  of  the  tilt-angle. 

4.2.3  Tool-Pin  Upper  and  Lower  Diameters.  The  effect 
of  the  tool-pin  size  (as  quantified  by  the  upper  and  lower 
diameter(s))  at  a  constant  value  of  the  tool-pin  taper  angle  of 
10°  can  be  inferred  by  comparing  the  results  displayed  in 
Fig.  14(f),  (a),  and  (g),  respectively.  The  tool-pin  upper 
diameters  in  these  three  cases  are:  0.012,  0.01,  and  0.008  m, 
respectively.  Examination  of  the  results  displayed  in  these 
figures  clearly  shows  that  the  highest  extent  of  marker-material 
mixing/dispersion  is  attained  in  the  case  of  the  largest  pin  size. 
Fig.  14(f)  and  that  this  extent  decreases  with  a  decrease  in  the 
pin  size.  This  finding  is  in  accordance  with  the  experimental 
observations  in  the  FSW  welding  practice  and  with  the  fact  that 
the  mechanical  energy  introduced  into  the  workpiece  through 
the  rotating  FSW  tool  increases  with  the  tool  pin  size. 

5.  Summary  and  Conclusions 

Based  on  the  results  presented  and  discussed  in  the  present 
work,  the  following  main  summary  remarks  and  conclusions 
can  be  made: 

1.  A  coupled  thermo-mechanical  Eulerian/Lagrangian  finite 
element  analysis  is  earned  out  to  investigate  material 


flow  and  mixing  during  the  Friction  Stir  Welding  (FSW) 
process.  Particular  attention  is  paid  to  the  effect  of  “weld 
pitch”  (defined  as  a  longitudinal  advancement  per  single 
revolution  of  the  tool),  toot  tilt-angle,  and  the  tool  pin 
size  on  the  flow  pattern  and  the  extent  of  material  mixing 
during  FSW. 

2.  The  results  obtained  show  that  FSW  does  not  cause  a 
large-scale  permanent  exchange  of  the  workpiece  material 
between  the  advancing  and  the  retreating  side,  i.e.,  major¬ 
ity  of  the  workpiece  material  within  the  weld  resides  on 
the  workpiece  side  (advancing/retreating)  on  which  it 
was  located  before  welding. 

3.  For  the  most  part,  FSW  causes  the  workpiece  material 
flow  in  the  horizontal  plane  (a  plane  parallel  with  the 
workpiece  backing  plate).  In  addition,  downward 
through-the-thickness  motion  of  the  material  and  the 
accompanying  mixing  are  found  to  be  quite  pronounced 
on  the  advancing  workpiece  side  while  this  component 
of  the  material  flow  was  not  very  prominent  on  the 
retreating  side. 

4.  The  desired  maximum  extent  of  material  mixing  has  been 
found,  within  the  explored  ranges  of  FSW  process 
parameters,  to  correspond  to  the  minimum  value  of  the 
weld  pitch,  an  optimum  value  of  the  tool  tilt-angle,  and  a 
maximum  value  of  the  tool  pin-size.  A  simple  rationale 
has  been  provided  for  these  observations. 
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